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Abstract 

We treat the vakonomic dynamics with general constraints within a new geometric 
framework which will be appropriate to study optimal control problems. We compare 
our formulation with Vershik-Gershkovich one in the case of linear constraints. We 
show how nonholonomic mechanics also admits a new geometrical description wich en- 
ables us to develop an algorithm of comparison between the solutions of both dynamics. 
Some examples illustrating the theory are treated. 
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1 Introduction 

As is well known, the application of tools from modern differential geometry in the fields 
of mechanics and control theory has meant a great advance in these research areas. For 
example, the study of the geometrical formulation of the nonholonomic equations of mo- 
tion has led to a better comprehension of locomotion generation, controllability, motion 
planning and trajectory tracking, raising new interesting questions in these subjects (see 
0. [|, rj, [H| [IT], [JT], and references therein). On the other hand, there is now 



a considerable amount of papers in which optimal control problems are addressed using 
geometric techniques. 



In this context, we present a unified geometrical formulation of the dynamics of nonholo- 
nomic and vakonomic systems. Both kind of systems have the same mathematical "ingredi- 
ents" : a lagrangian function and a set of nonintegrable constraints. But the way in which 
the equations of motion are derived differs. In the case of vakonomic systems, the dynamics 
is obtained through the application of constrained variational principle M. In particular, 
an optimal control problem can be seen as a vakonomic one. The term "vakonomic" ("vari- 
ational axiomatic kind") is inherited from Kozlov [|17|, who proposed this mechanics as an 
alternative set of equations of motion for a physical system under the presence of non- 
holonomic constraints. Nonholonomic equations of motion are deduced using d'Alambert 
principle when the constraints are linear or affine. 

The two approaches have deserved a lot of attention in recent years (see |1|, |7|, [19], |23|, 
|27| , |43j and references therein). Vakonomic mechanics (also called dynamical optimization 
subject to nonholonomic constraints) is used in mathematical economics (growth economic 
theory), sub-riemannian geometry, motion of microorganisms at low Reynolds number..., as 
we will illustrate in Subsections 3.1, 3.2 and 3.3, while, as we have mentioned above, non- 
holonomic mechanics has important applications to robotics, dynamics of wheeled vehicles, 
motion generation... 

Several authors have discussed the domains of validity of each mechanics [p], [H], p3f , and 
in fact, the question remains not completely closed. The solutions of these dynamics do 
not coincide, in general, though there are examples in which nonholonomic solutions can 



be seen as solutions of the constrained variational problem. In recent papers |9], ^3| the 
characterization of this situation has been studied. In [^j Lewis and Murray introduced 
the example of a ball on a rotating table and shown that the subset of solutions of the 
nonholonomic problem is not included in the set of vakonomic ones. In Favretti obtains 
conditions in some particular cases for the equivalence between the two formulations. 

Our project of unifying the comparison of both dynamics from a geometrical point of view 
has brought us to develop new geometric frameworks for vakonomic and nonholonomic 
mechanics, strongly inspired in the Skinner and Rusk formulation for singular lagrangians 
systems EDL By means of these approaches, we are able to compare them using an 



algorithm which gives rise, under appropriate conditions, to a final constraint submanifold 
describing the nonholonomic solutions which are also vakonomic. 

The paper is structured as follows. In Section 2, we obtain the equations of motion for 
vakonomic mechanics, assuming an admissibility condition, which permits us to present 
it using the restricted lagrangian to the constraint submanifold. Let us recall that from 
a geometrical point of view, the lagrangian is defined on the tangent bundle TQ of the 
configuration manifold Q, and M represents the submanifold of TQ determined by the 
vanishing of the nonholonomic constraints. We will treat here with arbitrary submanifolds, 
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that is, our constraints will be nonlinear in principle. It should be also remarked that we do 
not consider abnormal solutions. 

Section 3 is devoted to reformulate in geometric terms vakonomic mechanics. In this section 
we will use as ambient space the fibred manifold Wo = T*Q Xq M, which is in fact a 
subbundle of the Whitney sum T*Q @TQ (the phase space in Skinner and Rusk approach). 
Since T*Q is equipped with a canonical symplectic form we can induce a presymplectic 
structure u on T*Q Xq M. Moreover, we can consider the hamiltonian function Hy/ Q = 
i^i, ^2) — ^L-i where 71*1 and 7r 2 are the canonical projections, and L is the restriction of L 
to M. Then, we prove that the equations of motion of vakonomic mechanics are equivalent 
to solve the presymplectic hamiltonian equation 

Since the 2-form u is presymplectic, a constraint algorithm must be performed in order to 
obtain well-defined solutions of the dynamics. If the problem is consistent, we obtain a 
family of explicit solutions on the final constraint submanifold. In addition, a compatibility 
condition is found which characterizes when the first constraint submanifold W\ is symplec- 
tic, and therefore the algorithm stabilizes at the first step. Several applications to economy, 
locomotion and control theory, and subriemannian geometry are given in subsequent sub- 
sections. 

In Section 4, we compare our approach with that of Vershik-Gershkovich |E| for vakonomic 
systems with linear constraints. We prove that both are related by a convenient presym- 
plectomorphism, so that our approach could be consider as its generalization for the case of 
nonlinear constraints. 

Since we want to compare vakonomic and nonholonomic dynamics, it is necessary to con- 
struct a geometrical framework for nonholonomic mechanics using a close phase space. In- 
deed, in Section 5 it is proved that the nonholonomic dynamics lives on a submanifold M 
of Wq. In general, we have again a presymplectic system there and a constraint algorithm 
is needed to obtain the dynamics on the final constraint submanifold. 

In Section 6, assuming that the dynamics lives on W\ and M, respectively, we can compare 
their solutions by means of the mapping T : W\ — > M, (at,v) 1 — ► (Legi(v),v). Several 
illustrative examples are worked in order to show the different behaviours. It should be 
remarked that recent results by Favretti |§ are reobtained with weaker hypotheses. 
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2 Variational approach to constrained mechanics 



Let Q be an n-dimensional configuration manifold and L : TQ — > K. an autonomous 
lagrangian function. If (q A ), 1 < A < n, are coordinates on Q, we denote by (q A , q A ) 
the natural bundle coordinates on TQ in terms of which the tangent bundle projection 
t q : TQ — > Q reads as r Q (q A , q A ) = (q A ). 

Let us suppose that the system is subject to some constraints given by a (2n— Tridimensional 
submanifold M of TQ, locally defined by $ Q = 0, 1 < a < m, where $ a : TQ — > R. 

Along the paper, we will assume the admissibility condition for the submanifold M C TQ, 
that is, for all x G M, we have 

dim T X M° = dim S*T X M° , 

d 

where S = dq <S> ttt is the canonical vertical endomorphism (see [ETTI ) . This is equivalent 
oq A 

to say that the rank of the matrix 

<9($\...,$ m ) 
d(q\...,q n ) 

is m for any choice of coordinates (q A , q A ) in TQ. Consequently, by the implicit function 
theorem, we can locally express the constraints (reordering coordinates if necessary) as 

€ = * a (q A , q a ) , (1) 

where 1 < a < m, m + 1 < a < n and 1 < A < n. Then, (q A , q a ) are local coordinates for 
the submanifold M of TQ. 

We denote the set of twice differentiable curves connecting two given points x, y G Q as 

C 2 (x, y) = {c : [0, 1] — >Q\c is C 2 , c(0) = x and c(l) = . 
This set is a differentiable infinite-dimensional manifold 

Let c be a curve in C 2 (x,y). A variation of c is a curve c s in C 2 (x,y) such that Co = c. An 
infinitesimal variation of c is the tangent vector of a variation of c, that is, 

dc s {t) 



u(t) 



s=0 



The tangent space of C 2 {x, y) at c is then given by 

T c C 2 (x,y) = {u : [0,1] — >TQ / u isC 1 , G T c(t) g, «(0) = and «(1) = 0} . 
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Now, we introduce a special subset C 2 (x,y) oiC 2 (x,y) which consists of those curves which 
are in the constraint submanifold M 

C 2 (x,y) = {c G C 2 (x,y) | c(t) G M c(t) = M n r^^t)) , Vt G [0, 1]} . 

Let us consider the functional J defined by 

J : C 2 (x,y)-^R 

L(c(t))dt, 

Jo 

which we want to extremize among the curves satisfying the constraints imposed by M, 
ceC 2 (x,y). 

Definition 2.1 A curve c G C 2 (x,y) will be a solution of the vakonomic problem if c 

is a critical point of J\gi, xy \. 

Therefore, c is a solution of the vakonomic problem if and only if dj[c) • u — 0, for all 

u e T c C 2 {x,y). 



Remark 2.2 In this paper, we will assume that the solution curves c G C(x,y) admit non- 
trivial variations in C(x, y). These solutions are called normal in the literature, in opposition 
to the abnormal ones, which are pathological curves which do not admit nontrivial varia- 
tions PJ. Several investigators have shown the existence of C , stable under perturbations 



abnormal solutions 24, 29 



Now, we find a characterization for a curve c to be critical for the vakonomic problem. 

Proposition 2.3 A curve c G C 2 (x,y) is a normal solution of the vakonomic problem if 
and only if there exists /i : [0, 1] — > M. m such that 



d I dL\ dL 
dt \dq a J dq a " 

dL d^ 3 



I'c 



d fd^! a 



dt \ dq a J dq a 



(2) 



Q q a Vf3- dqa 

k q a = ^J a (q A ,q a ), 



where L : M —>■ R is the restriction of L to M. 
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Proof : The condition for a curve to be critical is 



= dj(c) ■ u = ^J(c s 



s=0 



dc s 

for any variation c s in C 2 (x, y) of c, where u = — - 

as 

Then, we have that 



s=0 



= — J(c s 
as 



In local coordinates, we obtain 



s=0 



£ L(c.{t)) dt 
Jo ds 



s=0 



l 



dL A dL . a dL d^ a a dL d^ c . „ . , 
u + 771 — u + — - A u + — — —u a ) dt 







dq A dq a dq A 



dq a dq a dq A dq a dq a 
dL dL d^ a 
dq a dq a dq a 



dL dL d^ a 
+ 



u A + 



u a dt 



(3) 



dL ^ A dL ^ a . hj 
o \ dq A dq a 



^From (|l|) we know that the infinitesimal variations u A , 1 < A < n, are not arbitrary. 
Consider the functions \i a defined as the solutions of the following system of first order 
differential equations 



dL 



Then, using the fact that u c 

d . 



d^ a 



dq a 

d^ a A . . 

ITa u + ' we S et 

dq A dq a 



1 < a < m . 











^ dq" 


V 


a dL d^ a 




dm* 


U dq* ' ^ <9g« 


U a + jJL a 


dq a 



or, equivalently, 



dL d 



u 



{ll a U C 



d^ c 



d^ a 



-U 



jjj a - 



-u 



dq a dt^ ' r dq a r dq a 
Substituting the last expression in (||) and integrating by parts, we obtain 



dj(c) • u 



dL 
dq a 



He 



'~d~t 



u a + 



dL 

dq a 



fJL c 



d^ a 



u a dt. 
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Now, since 



dL 



dq a 



' dq a 



u 



dL 

dq a 



' dq a 



using again integration by parts, we can write 

-i 







o 



u 




jjj a ■ 



d^ c 
dq a 



dL a d^ a 
dq a ^ dq a dt \ dq a 



d ( dL d^! a 



dq a 



u a dt . 



u 



As the infinitesimal variations u a are arbitrary, the fundamental lemma of the Calculus of 
Variations applies and we can assert that dj{c) ■ u — if and only if c and fi a satisfy 
equations (0). \qed\ 



Remark 2.4 The usual way in which the equations of motion for vakonomic mechanics are 
presented is the following 



d f dL 



dL 



, d$ a 

Art „ . * + A c 



dt \dq A J dq A dq A 
$> a (q,q) =0, 1 < a < m , 



d fd$ a 



dt \dq A 



d$ c 
dq A 



(4) 



where $ Q = ^ a -q a and A c 



dL 

dq a 



— Ha, 1 < a < m . Observe that, in contrast to equations 



(g), equations (|J) are expressed in terms of the ambient lagrangian L : TQ — > R. Equations 
(0) stress how the information given by L outside M is irrelevant to obtain the vakonomic 
equations, contrary to what happens in nonholonomic mechanics (see Section || below). 

Equations (f|) can be seen as the Euler-Lagrange equations for the extended lagrangian 
C = L + A a $ Q . We will not follow this approach here, which has been exploited fruitfully 
in 0, [TBI |2"%fl. Finally, note that if we consider the extended lagrangian Aq£ + A a $ a , with 



A = or 1, then we recover all the solutions, both the normal and the abnormal ones 



3 Geometric approach to vakonomic mechanics 



We will develop a geometric characterization of vakonomic mechanics following an approach 
similar to the formulation given by Skinner and Rusk II for singular lagrangians (see 
also |§ [12], ^5|). This characterization is specially interesting, since enables us to study 
both linear and nonlinear constraints in an intrinsic way. We will show in the examples the 
utility of this formulation. 
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Consider the Whitney sum of T*Q and TQ, T*Q © TQ, and its canonical projections 

pri : T*Q © TQ — > T*Q , pr 2 : T*Q © TQ — ► TQ . 

Let us take the submanifold Wq = pr^" 1 (M), where M is the constraint submanifold, locally 
determined by the constraint equations $° = 0, 1 < a < m. We will denote Wq = T*QxqM 
and 7Ti = pri^, 7r 2 = p^i^. Now, define on T*Q Xq M the presymplectic 2-form uj = 
7t*o;q, where cjq is the canonical symplectic form on T*Q. Observe that the rank of this 
presymplectic form is equal to 2n everywhere. Define also the function 

H\v Q = (7ri,7r 2 ) - , 
where L : M — > R is the restriction of L to M. 

If are local coordinates on a neighborhood U of Q, Q ,a ) coordinates on TUOM and 
(<? A ;Pa) the induced coordinates on T*U, then we have induced coordinates (q A ,PA,Q a ) on 
T*U Xq (TU H M). Locally, the hamiltonian function ifyK reads as 

^Wo(/, PA, (f) = Paf + Pa ^ a ~ L(q A , q a ) , 
and the 2-form uj is u = dq A A dp a. 

Now, we will see how the dynamics of the vakonomic system (0) is determined by studying 
the solutions of the equation 

i x uj = dH Wo . (5) 
Thus, we are justified to employ the following terminology: 

Definition 3.1 The presymplectic hamiltonian system (T*Q Xq M,uj,Hw ) will be called 
vakonomic hamiltonian system. 

Being the system (T*Q Xq M,uj, H Wo ) presymplectic, we apply to it the Gotay-Nester's 
constraint algorithm [[iOj [TT| . First we consider the points W\ of T*Q Xq M where (y) has 
a solution. This first constraint submanifold is determined by 

W\ = {iG T*Q x Q M I dH Wo (x)(V) = 0, e kerc;(a;)} . 
d 

Locally, keru = (— — ). Therefore, the constraint submanifold W\ is locally characterized 
oq a 

by the vanishing of the constraints 

d^ a dL 

Va=Pa+Pa W - Wa =0, 
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or, equivalently, 



dL d^ a 
Va = 77T7 - VolTttt , m + 1 < a < n . 



(6) 



dq a dq a 

Expanding the expressions in equation (||) and equating coefficients, we obtain that the 
equations of motion along W\ are 

. A dH Wo 



Pa 



dp A 
dH Wo 

dq A 



which is equivalent to 



Q 

Pc 



8L d^ p 
Pfs- 



d dL d^ c 
dt \ W~ Pa ~dq~ a 



dq a 
dL 



Pf3- 



dq a 



(7) 



dq a dq a 



Observe that these equations are precisely the vakonomic equations of motion (gj), where 
now p a = fi a , 1 < a < m. 

Remark 3.2 The momenta p a , 1 < a < m, play the role of the Lagrange multipliers, but 
they do not have any physical meaning (see |f42|). 



Therefore, a vector field X solution of equation (||), will be of the form 



X 



„ , d I d 2 L 
9 I— + 



dq a ' \ dq a dq b 



Pi 



d 2 W \ d 



+ Y a [ — + 



dq a dq b J dp b 
d 2 W \ d 



d_ 

dq c 



+ 



d 2 L 



dq a dq b Vl 'dq a dq b J dp b 
( d d^ a d 



d 2 W \ d 



dL 

dq a ' \dq a dq b ^ 7 dq a dq b I dp b I \ dq a 13 dq a ) \dp a dq b dp b 



Nevertheless, the solutions on W\ may not be tangent to W\. In such a case, we have 
to restrict W\ to the submanifold W2 where these solutions are tangent to W\. Proceeding 
further, we obtain a sequence of submanifolds (we are assuming that all the subsets obtained 
are submanifolds) 



W, 



W 2 ^W 1 ^W = T*Q x Q M . 



Algebraically, these constraint submanifolds may be described as 

Wi = {xE T*Q x Q M I dH Wo {x){v) = , Vv e T X W^ } , i>l 
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T*Q © TQ 




Geometric formulation of vakonomic mechanics 



where 

T X W^ = {ve T X (T*Q x Q M) | u(x)(u,v) = , Vw G T x W t ^ } . 

If this constraint algorithm stabilizes, i.e., if there exists a positive integer fceN such that 
Wk+i = Wk and dimM 7 *; 7^ 0, then we will have obtained a final constraint submanifold 
Wf = Wk on which a vector field X exists such that 

(i x uJ = dH Wo )^ Wf . 

Note that on Wf we will have an explicit solution of the vakonomic dynamics. A very 
important particular case is when the final constraint submanifold is the first one, Wf = W\. 
Observe that the dimension of W\ is even, dim W\ = In. In the sequel, we will investigate 
when this constraint submanifold is equipped with a symplectic 2-form in order to determine 
an unique solution of the vakonomic equations. Obviously, this geometrical study is related 
with the explicit or implicit character of the second order differential equations obtained in 

(!)• 

Denote by the restriction of the presymplectic 2-form to to W\. 



Proposition 3.3 (Wi,uwi) is a symplectic manifold if and only if, for any choice of coor- 
dinates (q A ,pA,q a ) on T*Q x Q M, 

( d 2 L _ d 2 ^ a . 

for all point in W\. 
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Proof: luw i is symplectic if and only if 

T x W 1 n(T x W 1 ) ± = 0. 



(8) 







for all x G W\. Condition (|8|) is satisfied if and only if the matrix dip a (— — ) is regular, that 



is. 



d 2 L d 2 ^ a 
dq a dq b ^ a dq a dq b ' 



for all x G W\. 



\QED\ 



In this case, the equations of motion (|?p are rewritten as the following system of algebraic 
and explicit differential equations 



q" 



Po 

q a 



ty Q {q A ,q a ) 



dL 

dq a 

-C ab 
dL 



Pp 



dq l 



where 



+ Pc 



C, 



dm? 

d 2 L 
QqAQgb 



- Q A Pa 



dL 



d 2 ^ a 



QqAQgb 



(9) 



dq b \ dqi 

d 2 i 



Pf3 



dq'y ) dq b 
d 2 ^ a 



ab 



QgaQqb 

and (C ab ) denotes the inverse matrix of (C a &)- 



-Pa 



dq a dq 



(10) 



Remark 3.4 The characterization found in Proposition |3.3| for the symplecticness of the 
manifold (Wi,ujwi) implies that the constraint equations 

d^ a dL 



<fia = Pa +Pa- 



—- = , m + 1 < a < n . 



dq a dq a 

define locally the variables q a , m + 1 < a < n, by the implicit function theorem. That is, 
we have 

q a = S a (q A ,p A ) , m + 1 < a < n . 

Therefore, we can also consider local coordinates (q A ,PA) on W\. In such a case, the sym- 
plectic form has the following local expression 

uj Wi = dq A A dp a , 

and the restriction of the hamiltonian Hy/ Q to W\ 

H Wl = PaS a + p a ^ a - L(q A ,PA) , 
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where L(q A ,pA) = L(q A ,s a (q A ,pA))- Consequently, equations @ are rewritten in hamilto- 
nian form as 



Pa 



% Vl ("J 



dq A 

This choice of coordinates is very common in optimal control theory. 

Now, observe that, if the constraints are linear on the velocities, we can write 

q a = ^ a (q)q a . 



Then, from Proposition 3.3, u^i is symplectic if and only if 



Proposition 3.5 Suppose that the constraints are given by 

q a = ^ a a (q)q a , 1 < « < m, 

and the lagrangian L is regular. Denote by (W AB ) the inverse matrix of the Hessian matrix 
of L. In this case, u>w 1 is symplectic on W\ if and only if the constraints are compatible, 
that is, the matrix whose entries are 

is nonsingular. 



Proof : See the geometrical proof of Theorem IV. 3 in reference \qed\ 



Remark 3.6 The compatibility condition guarantees the existence and unicity of the so- 
lutions for the nonholonomic problem with lagrangian L and constraint submanifold M 



ia 39 . 



3.1 Applications to economy 

The variational calculus is an indispensable tool in many classical and recent economic 



papers [|13|, |(| |35|, |37], |38[. In fact, a typical optimization problem in modern economics 



deals with the problem of maximizing or minimizing the functional 



f 

Jo 



D(t)U[f(t,k,k)} dt 
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subject or not to constraints. Here, D(t) is a discount rate factor, U an utility function, / a 
consumption function and k the capital labor-ratio. It is usual to find dynamical economic 
models with nonholonomic constraints. For instance, the revision of the expected rate of 
inflaction may be expressed in terms of the nonholonomic constraint 

7T =j(p- n ), 0< j < 1 

where 7r and p are the expected and actual rates of inflaction, respectively. 

In economics, it is also very common to deal with an explicit dependence of the time. We 
plan to extend the geometric formulation of vakonomic dynamics to the non-autonomous 
case in a forthcoming paper. 



Example 3.7 (Closed von Neumann System ||37| , |38|| ) Consider the transformation func- 
tion which relates n capital goods K\, K 2 , . . . , K n and the net capital formations K\, K 2 , . . . , K n 



F(K h ...,K n ,k x ,...,k n ) = K^K? ...K^-\kl + ... + k 2 n 
with a,\ + «2 + • • • + On = 1- The von Neumann problem is to maximize 

K n dt subject to F(Ki, . . . , K n , K\, . . . , K n ) = 



1 1/2 





and appropiate initial conditions. 

Applying our formulation it is possible to write this problem as a presymplectic system on 
Wo = -R 3n_1 . The constraint F = can be rewritten as 

\ 1/2 

'2ai 



■^(K 1 ,...,K n ,K 2 ,...,K n ). 



K 1= \ Kr,...,K^-J2^ 

i=2 

Taking coordinates (Ki, . . . , K n , K 2 , . . . , K n , P 1 , . . . , P n ) we have that 

n 

dKi A dP i 



UJ 



1/2 



h Wo = pt ^ + pl -[ K i ai K 2 a2 ■ ■ ■ K an - J2 - K * 



i=2 



i=2 



Applying the Gotay and Nester algorithm new constraints arise 

-1/2 



P i = P l Ki Kl ai K 2 2 a2 ■ ■ ■ K 2 n an - K , 

\ i=2 

P n = 1 + p l K n ( K 2ai K 2a2 ■ ■ ■ K 2 n an - K , 



2<i<n-l 



\ 1 ^ 

i 

i=2 
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Therefore, from (^) the initial system is determined solving the following n differential 
equations on the variables (Ki, . . . , K n , K 2 , . . . , K n , P 1 ) 



= P x k;G 



where 



+P 1 



K**).G 



2ai-l 



G(Ki, . . . , K n , K 2 , . . . , -ft'n 



(12) 



K*-)) G + k4(G) 



, 2 < % < n 



V(K h ...,K n ,K 2 ,...,K n 



3.2 Principal kinematic locomotion systems 

This kind of systems includes the motion of inchworms, paramecia, mobile vehicles, robotic 
snakes, etc [[14], [T^, ^3j. The study of the motion relies on the simple fact that the process 
of locomotion can be divided into internal (shape) variables and position (group) variables. 
The internal variables are assumed to be directly controlled. The motion of these variables 
couple to produce a net change in the position and orientation of the moving body. Moreover, 
locomotion systems are characterized by the fact that the constraints are usually invariant 
with respect to the group action. Consequently, one is provided with a useful mathematical 
structure to work with: a principal connection 7 on a principal fibre bundle, 7r : Q — ► Q/G. 
The constraint submanifold is precisely the horizontal distribution 7i of 7. 

In the following, we will investigate optimal control problems for these systems. Take local 
coordinates in Q, (r,g), where r are coordinates in the base manifold Q/G and g in the Lie 
group G. Let us assume that we are given a quadratic cost function C, locally expressed as 

C(r, r) = -C a b(r)r a r b and define the functional 
2 



j 



C(r,f)dt. 



The cost function C depends only on the shape variables, which corresponds to calculating 
the cost of the control effort. Then, the optimal control problem is to obtain the inputs that 
will minimize J, while steering the state from (ro,<7o) to (ri,gi). 

^From the vakonomic point of view, we have the lagrangian L = C, the constraint sub- 
manifold M = 7i and the infinite-dimensional manifold C 2 ((r , g ), (r 1; g\)). We will not 
go into depth here with the mathematical structure associated to the connection 7 and the 
symmetries G. This will be a subject of future research. 
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Example 3.8 (Locomotion at low Reynolds' number |fj"5| , |31|| ) The kinematic connec- 
tion for the paramecia can be determined by examining the Stokesian flow around a de- 
formable cylindrical body [15 1 . We parametrize the body in polar coordinates (R, 9) as 



R 



e(ki(t) cos 20 + k 2 (t) cos 39) 



The shape variables for this system are r = (hi, k 2 ). Let x denote the motion of the centroid 
of the body in the direction given by 9 = 0. Symmetry arguments show that all resultant 
motion must be directed along this ray. Therefore, we have a principal bundle 7r : R 3 — > M 2 
with structure group G = R. In ||15|| , it is shown that the viscous connection can be written 
to first order as 

e 2 • 

7 = x + — {k 2 k\ + 2k x k 2 ) . 

We consider the optimal control problem associated with the simple quadratic cost function 
C = /c 2 ~\~ A^2* 

Taking coordinates (ki, k 2 ,x, ki, k 2 ,pi,p 2 ,p 3 ) in W = M Xq T*Q, we have that 



H 



uj = dki A dpi + dk 2 A dp 2 + dx A dps 

e 2 ... 

Wo = Pih +P2h -P3 — (k 2 ki + 2k ± k 2 ) - k\ - k\ 



Applying the Gotay-Nester algorithm, new constraints arise 

e 2 

Pi = -^k 2 p 3 + 2ki , 
e 2 

P2 = -7^ k lP3 + 2k 2 ■ 



Now, equations (0) read as 



x 

P3 

h 
h 



-j(k 2 k 1 + 2k 1 k 2 ) 



P317h 



6 I 
-P3Trh 



(13) 



These equations can be easily integrated. Setting a = p 3 —, we have 



^4 

ki(t) = B cos at + C sin at 

a 

k 2 (t) = B sin at — C cos at + D 



for some constants A, B, C, D, which is the same result obtained in [31 
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3.3 Vakonomic mechanics and sub-riemannian geometry 



Let Q be an n-dimensional manifold with a smooth distribution V of constant rank n — m. 
A sub-riemannian metric |j| on T> is a smoothly varying in q positive definite quadratic 
form g q on V q . A piecewise smooth curve 7 in Q is called admissible if j(t) G V. We define 
the length of such a curve in the usual way 



Z( 7 )= / y/gW)Mt))dt. 



The sub-riemannian distance between two points x, y G Q is defined as 

d(ar,j/) = inf 7 (/(7)) , 

for all admissible curves 7 connecting x and y. The distance is taken infinite if there is no 
such a path. 

A curve which realizes the distance between two points is called a minimizing geodesic. 
It is easy to show that 7 is a minimizing geodesic if it minimizes the functional 

among all the admissible curves with the same endpoints. 

Let fix, . . . , fi m be a basis of 1-forms for the annihilator T>° . Then an admissible path must 
verify the nonholonomic constraints 

^(7) =0, 1 < i < m . (14) 

Thus, we see that the problem of finding minimizing geodesies corresponds exactly to the 
problem of solving the vakonomic problem determined by the restricted lagrangian L = ^g 
and the nonholonomic constraints (14). Note that, as the constraints are linear and g is 
positive definite, the Gotay-Nester's algorithm always ends in the first step, Wf = W\. 

Example 3.9 (Sub-riemannian geometry: The Martinet case [H) Let U be an open 
set of 1R 3 containing and T> the distribution on U determined by annihilation of the Mar- 
tinet 1-form 

y 2 

fi = dz — 2 

y 2 

which determines the constraint function z = — x. Moreover, consider the restricted la- 

2 

grangian 

L(x,y,z,x,y) = -(x 2 + y 2 ) . 
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This corresponds to the flat metric case in |J. 

The normal minimizing geodesies for this problem are determined by solving the presym- 
plectic hamiltonian system 

i x uj = dH Wo , 

where, locally, 

uj = dx A dp x + dy A dp y + dz A dp z , 
H Wo = xp x + yp y + ^-xp z -^(x 2 + y 2 ) . 

Applying the Gotay-Nester constraint algorithm we obtain the new constraints 

. y 2 

Px = X - —p z , 

Py = V 

In this particular case, the equations of motion ([?]), in coordinates (x,y, z,x,y,p z ), are 

y 2 . 

z = -x 

Pz = 

d (. y 2 \ ( 15 ) 



i x p z = 

r// V 2 y J 

y = -yxp z 

which are obviously integrable by quadratures (compare with 0). 

4 Vershik-Gershkovich and vakonomic hamiltonian ap- 
proaches compared 

In the precedent section, we have found an intrinsic geometric approach to vakonomic dy- 
namics. It is possible to give an alternative geometric formulation of the vakonomic equations 
of motion related to the one of Vershik and Gershkovich p31 . A key element to obtain this 
alternative description will be the next fibred diffeomorphism 

F : T*Q © TQ — ► T*Q@TQ 

(a,v) i — ► (a - Leg L (v),v) , 



for any a G T*Q, v G T X Q and x G Q. Here, Leg L : TQ — > T*Q denotes the Leg 

>n 'c 
dl 

dq 



endre transformation associated to the lagrangian L, which in local coordinates reads as 
dL 

Leg L {q A ,q A ) = {q A , -^tj). It is clear that F(T*Q x Q M) = T*Q x Q M. We will see how in 
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the case of linear constraints, we "recover" the Vershik-Gershkovich formulation. As a by- 
product, we will have obtained a generalization of their formulation to the case of nonlinear 
constraints. 

Consider on T*Q © TQ the presymplectic 2-form Q = pr*ujQ. Let ujl = —dS*dL be the 
Poincare-Cartan 2-form on TQ associated to L : TQ — > K and El its energy function. Take 
also the presymplectic 2-form pr\iOL on T*Q © TQ, and define the functions 

H = (pr 1 ,pr 2 ) -pr 2 L, 
H = (pr 1 ,pr 2 ) - pr* 2 E L • 



Lemma 4.1 The diffeomorphism F : T*Q © TQ — ► T*Q © TQ zs a presymplectomorphism 
from (T*Q ®TQ,tt) onto (T*Q ® TQ,tt + pr* 2 tu L ) , i.e., F*(tt + pr*co L ) = ft. Moreover, it 
verifies that F*H = H. 



Proof : F is clearly invertible with inverse 

T- 1 : T*Q © TQ — ► T*Q © TQ 

(a,v) i — ► (a + Leg(v), v) . 

A direct computation shows that H o T _1 = H. Moreover, in local coordinates, 

dL 



(F- l )*(dq A A dp A ) = dq A A 



dL 



dq A dp a + dq A d 



dq J 



a y ' 



which implies F*(Q + prluj L ) = fi. |Qgp| 

Denote by j : T*Q XqM <^-> T*Q@TQ and i : M ^ TQ the respective canonical inclusions. 
Let us define u; = + pr^L)- Since pr 2 o j = i o 7r 2 , we have that 



W = W+ JO 7T 2 ) C^L . 



Proposition 4.2 T/ie solutions of the equations 

i x uj = dH Wo , (16) 

and 

l Y Uj = d(fH) , (17) 

are F\w -related, that is, if x G T*Q Xq M is a point where there exists a solution Y of 
equation fcF\ ) then TT -1 (Y) is a solution of equation fllb] ) at F~ l (x) and, conversely, if X 
is a solution of equation /\ldj) at T _1 (x) then TF(X) is a solution of equation fcDQ at x. 
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Proof : It readily follows from Lemma 4.1. \qed\ 



An inmediate consequence is the following 

Corollary 4.3 F preserves the constraint submanifolds provided by the presymplectic sys- 
tems (T*Q Xq M,uw , Hw ) and (T*Q Xq M,u),j*H). That is, if 

...^W k ...^W 1 ^W = T*Q x Q M and 

...^P k ...^P 1 ^P = T*Qx Q M, 

are the sequences of submanifolds generated by the Gotay and Nester's algorithm for the first 
and the second presymplectic hamiltonian system, respectively, then 

F t = F m : Wi _> i* , 

are diffeomorphisms for all i . 



In conclusion, Proposition ^]2| and Corollary |4.3| show that it is equivalent to solve the vako- 
nomic hamiltonian equations (fL6|) as in Section | or equations (|17|). Locally, if (q A (t), pA(t), q a {t)) 
is an integral curve of X then 

(q A (t),p A -i^{q B (t),q b (t)),q a (t)) 

is an integral curve of Y . 

Next, we will study solutions of equations (|17]) from a local point of view. First, notice that 
it is clearly equivalent to solve equations (|T7|) or the following system of equations on the 
ambient space 

(i Y {n+pr*u; L ) = dH)^ QxQM 

^\t*Qx q M ^ T(T*Q Xq M) . 

Now, if (q A , X A , q A ) are local coordinates on T*Q®TQ and (q A (t), X A (t), q A (t)) is an integral 
curve of Y, then, since 

and 

Cl + pr\u) L = dq A A d\ A + dq A A d 

we have that 



dL 
dq A 



d 2 L d 2 L dL 
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When restricting this formula to T*Q Xq M, we must use that q a = ^ a (q A ,q a ) and then, 
d^ a A d^ a 

dq a = . dq + dq a . Therefore, the former system is written as 



dq 1 



dq a 



d ( dL 



dL 



dt \dq A J dq 



dq 1 



-\Adq J 



\M 



^ j A 

K-^-rdq 
dq A 



Xa + A, 



d^ a 
'~dq^ 



dq a = 



(18) 



Moreover we have that 



Xa ~Jt 



-A, 



d^ a 
~dq^ 



"Ar 



dq a 



d fd^> c 

Ary — 



dt V dq a 



After some computations, the system becomes 



d f dL 



dt \dq 



dL 
dq 1 



Art „ . . h A a 



a dq J 



d 


(d$ a \ 


d$ a ~ 


dt 


ydq A ) 


dq A _ 



$ a (q,q) =0, 1 < a < m, 

where now g) = ^/ a (q A ,q a ) — q a , 1 < a < m. These last equations are the clas- 

sical equations of motion for a vakonomic system or for the dynamic optimization under 
nonholonomic constraints (see equations (HI)). 



4.1 Vershik-Gershkovich approach 



In |^3[ , Vershik and Gershkovich gave a formulation for the "nonholonomic variational prob- 
lem" , i.e., the vakonomic problem, within the framework of the so-called mixed bundle, which 
we briefly review in the following. 

If T> : Q — > TQ is a differentiate distribution along Q then, the mixed bundle over Q 
associated to T> is given by D © T>°, where T>° is the codistribution annihilating T>. This is, 
the fibres of V © V° — ► Q are V q © V° q . 

Let {& a (q A ,q A ) = ^"(q)q a — q a , 1 < a < m } be a set of independent functions whose 
annihilation defines the distribution T> and let {i] a = ^/^dq a — dq a , 1 < a < m } be the 
corresponding basis of T>°. Regarding T> C TQ as the set of admissible velocities, Vershik 
and Gershkovich write the equations of motion @ for the vakonomic problem (L, T>) as 
follows 

' {i{w)-W} iqA ^ Kt+K[Hith (19) 

< q, r] a >= , 1 < a < m . 
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In this particular case, we obtain that Pi, the first constraint submanifold for the presym- 
plectic hamiltonian system (T*Q Xg M,u>, j*H), is just T>° © V, since 

A a + KK = 0, 1 < a < m , 

from equations (|18|). 

If (Pi = V° © V, oupj) is a symplectic manifold (see Proposition |Q|) , then the equations of 



motion (p79f ) determine a unique vector field on D © Z>° and the Lagrange multipliers X a are 
coordinates in T>° with respect to the basis r] a . 

Consequently, the geometrical picture we have developed in Section [3] is equivalent to 
Vershik-Gershkovich approach. As said above, we have obtained a generalization of Vershik- 
Gershkovich formulation to the case of nonlinear constraints, just "translating" things from 
our approach by the diffeomorphism F. 

In the nonlinear case, under the admissibility condition, one can verify that the first con- 
straint submanifold Pi = F(Wi) can be identified with the manifold S*{TM°) Xq M. In 
fact, we have that S*{TM°) is generated by the 1-forms 

S*d$ a = dq a - ^—dq a , 1 < a < m . 
oq a 

If (\A,q A ,q a ) G Pi, then the 1-form \Adq A is a linear combination of the 1-forms S*d& a in 
the following manner 

\ A dq A = \ a S*d$ a . 



5 Geometric approach to nonholonomic mechanics 

A nonholonomic lagrangian system consists of a lagrangian L : TQ — > IR subject to non- 
holonomic constraints defined by m local functions Q a (q A ,q A ), 1 < a < m. The equations 
of motion for nonholonomic mechanics are derived assuming that the constraints satisfy 
d'Alembert's principle, in the linear or affine case. In the nonlinear case, it does not exist 
an unanimous consensus about the principle to adopt |^7|, . The most widely used model 
is the Chetaev's principle and it will be assumed in this paper. The equations of motion are 
then given by 

d ( dL\ dL d<S> a 

A Ql — , (20) 



dt \ dq A J dq A a dq A 

together with the algebraic equations $ a (q A ,q A ) = 0. The functions X a , 1 < a < m, are 
some Lagrange multipliers to be determined. 
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As in the vakonomic case, we assume the admissibility condition, so it is possible to write 
the constraints as q a = ^ a (q A , q a ) , where 1 < a < m, m + 1 < a < n and 1 < A < n. 

The study of nonholonomic systems in the realm of Geometric Mechanics has been an 
active area of research in the last years (see, for instance, [19[ and references therein). The 
nonholonomic equations of motion can be written geometrically as 

(i r u L -dE L ) lM ES*(TM°), 
T ]M eTM, 

where the subbundle S* (TM°) of T*TQ along M represents the constraint forces. 

Nonholonomic mechanics also admits a nice geometrical description on the space T*Q ®TQ 
inspired in the one by Skinner and Rusk |40], |4lfl . In addition, this description will be ap- 
propiate to compare the solutions of the dynamics between the vakonomic and nonholonomic 
mechanics. In the following, we will prove that equations fl2~l~|) are equivalent to the next 
ones 

(i x tt - dH) lT » QxgM G F° , 
X\t*qx q m £ T(T*Q Xq M) , 
where Q is the presymplectic 2-form Q = pr\ujQ on T*Q © TQ, H the hamiltonian function 
H = (pr u pr 2 ) -pr*L and F° the subbundle of T*(T*Q © TQ) along T*Q x Q M defined by 
F° = pr*(S*(TM°)). 

Indeed we have in local coordinates 

= dq A A dp A , 

itt • a i i ■ a A dL , A 

dH = q dp A + p A dq - Tr^dq - —jdq , 

oq A oq A 

and F° is generated by the 1-forms 

l^A d< i = l^r^ ~dq a ,l<a<m. 
oq A oq a 

d d d 

If X = X A — —r + Y A — — j + Za~^ — was a solution of equations then we would have 
oq A oq A opA ~ — 

X A = q\ 

dL d$ a 

Za - W + a W ( } 



along with the constraints 



dL 

dqf 

§ a (q A 7 q A ) = 0. (24) 



Pa-tta = 
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Observe that these constraints determine the submanifold M of T*QxqM. The submanifold 
M is diffeomorphic to M since 

M — > M 

m i — > (LegL{m),m) , 

is a diffeomorphism. M is the first constraint submanifold provided by the constraint algo- 
rithm applied to equations fl2"2]). This algorithm will lead to a final constraint submanifold 
on which there exists a well-defined dynamics, at least in case the given problem is consistent 
(see 



19|1). Obviously, equations fl23|) and ( f24|) are equivalent to the nonholonomic equations 
of motion (f20|). 

In terms of the \I/ Q 's the above equations are written as 

X A = q\ 

dL d^ a 

dq a a dq a ' 



dL 



A 



along with the constraints 



Pa ~ 



Therefore, a solution X of 

d 



X 



dL 
dq A 

q a -V a {q A ,q a ) 
2|) is of the form 

d 2 L 



Q 

+Y a 



d^ a d 
dq a dq a dq a 
d d^ a d 



+ 



0. 
. 

d^ a d 2 L 



(25) 







dtp 
d 



+ 



dq"< dq c 



+ 



+ 



d 2 L d^ a d 2 L 



dq A dq^ 
d 2 L 



+ 



d 



+ 



Qqi dq A dq a ) dp a 
d^ a d 2 L \ d 



dq a ' dq a dq a ' \dq A dq a dq a dq A dq c 

Under the regularity assumption, which means that the matrix 

d 2 L .J dL\ d 2 ^ a 



dp A 



C, 



ab 



QqaQqb 



— I 



dq a ) dq a dq b ' 



(26) 



is invertible (see |39|]), there is an unique solution of the dynamics on M. In particular, after 



some computations, we obtain 

yd Q ab 

_dL_ 

dq 1 



+ i 



d 2 L 


-4 A i*{ 


QqAQqb 


■( dL ) 


fd^ a 


\dq a J 


V dq b 



dL\ d 2 ^ a 



dq a I dq A dq 



d^ c 



;b 



dq 



(27) 
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where i : M — > TQ is the canonical inclusion and C a& the inverse matrix of C^. 

Taking coordinates (q A , q a ) on M, the equations of motion for a nonholonomic system will 
be 



^ a (q A ,q a 
_ C ab 

dl 



q 



dq 1 



+ i 



d 2 L 




QqAQgb 






\dq a J 


V dq b 



dq a J dq A dq b 



(28) 



dq b 



Compare them with equations 



6 Vakonomic and nonholonomic mechanics: Equiva- 
lence of dynamics 

In this section, we shall investigate the relation between vakonomic and nonholonomic dy- 
namics. 

Consider a physical system with lagrangian L : TQ — > M and constraint submanifold M C 
TQ. Let us assume that the vakonomic problem lives in the first constraint submanifold, 
Wi, and that the nonholonomic one lives in M (this will be the case if the constraints are 
linear and the admissibility and compatibility conditions are satisfied). As a consequence, 
we have well defined vector fields X V k on W\ and X n h on M . 

It is clear that the mapping (vt 2 )|vki : W\ — > M is a surjective submersion and that we can 
define the mapping T : W\ — > M as 

T : Wi — > M 

(a,v) i — ► (Leg L (v),v) 

In coordinates, T reads as T(q A ,q a ,p a ) = (q A ,q a ). 

Our aim is to know whether, given a nonholonomic solution, we can find initial conditions in 
the vakonomic Lagrange multipliers, p a , so that the curve can also be seen as a vakonomic 
solution. In order to capture the common solutions to both problems, we have developed 
the following algorithm. It is inspired in the idea of the T-relation of X v ^ and X n h and the 



constraint algorithm developed by O. Krupkova [18|]. If both fields were T-related, then the 
projection to M of all the vakonomic solutions would be nonholonomic. So, selecting the 
points in which they are related, we are picking up all the possible good candidates. We 
write W\ = So and define 



S, = {w G S | T w T(X vk (w)) = X nh (T(w))} . 
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In general Si is not a submanifold. If Si = 0, there is no relation between the vakonomic 
and nonholonomic dynamics. 

If Si 7^ 0, we apply the following algorithm: 

• Step 1: For any w G Si, consider C( TO ) = UjC^j, the union of all connected submani- 
folds C(u,)i of maximal dimension lying in Si, contained in a neighbourhood U of w and 
passing through w (maximal dimension means that if N is a connected submanifold 
lying in Si fl U passing through w and C( w yi C A" then C( TO )j = AT). 

Suppose that Cr w ) 7^ {^}- For each i we consider the subset of Ci w )i 

C{ w )i = {v E C( w )j | X vk {v) e T v C( w )i} . 

If Cr w )i = Ct w )i then we call the submanifold C( w )i a final constraint submanifold 

at w. If C( w )i = 0, we exclude C( w )i from the bunch C^ w \. If C C^ w yi C C( w )j, then 
we proceed to the next step. 

• Step 2: Repeat the Step 1 with instead of Si. 

After sufficient steps of this algorithm we obtain a bunch of final constraint submanifolds at 
w or we find that there is no final constraint submanifold passing through w. Collecting all 
the points where there exists a bunch of final constraint submanifolds we obtain the subset 
where there is equivalence between vakonomic and noholonomic dynamics. 

Suppose that the constraints $ Q! , 1 < a < m, are linear on the velocities so we can write 
them as 

In such a case, the matrices C and C defined in (pi]) and (p6|), respectively, are the same 
(even for constraints affine on the velocities). 

Proposition 6.1 5*1 is locally chararacterized by the vanishing of the n — m constraints 
functions on Wi 

■a ( -*9L \ 

9b = q \ Pa - l w) 

Proof : The comparison between the vector fields X vk and X n h consists of taking the 
difference between q a, s in the expressions (^) and ( pgp and equating the result to zero. 

\qed\ 



dq a dq b 



a dqP 



h dq? 



m + 1 < b < n 



(29) 
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Consider the local projection p(q a , q a ) = (q a ) and the connection T on p such that the 
horizontal distribution 7i is given by prescribing its annihilator to be 

H° = (dq a - ^ a a dq a , l<a<m) . 

Then the curvature R of this connection (see ||21|| ) is given by 



d d d 



dq a) dq b ' ao dq a ' 
where 

R\ = b - a + b - a 

ab g q a g q b a b Q q p • 

We say that V is flat if the curvature R vanishes identically. The tensor R measures the lack 
of integrability of the horizontal distribution 7i, which in our case is the constraint manifold. 

Then, we can write the constraints determining Si as 

9b = q a (p a ~ **|^J Kb, m + l<b<n. 

^From this expression we obtain that if the constraints are holonomic, then R = and the 
final constraint submanifold is equal to So = W\. Therefore, every nonholonomic solution is 
also a vakonomic solution. Indeed, equations (0) will read as 

d ^ a -a 

i>» '"■Ihr' 1 " ■ (30) 

Q q a a Q qa 

The first and the third set of equations determine the trajectory in M. The Lagrange 
multipliers p a are determined by the second set of equations once we know the solution in 
M. This is the typical behavior of the holonomic case f23|. 

But, in general, for linear constraints, the first constraint subset in the algorithm is deter- 
mined by 

Si = {9b = 0, m + 1 < b < n} , 
dL 

where gb(q A , q a , Pa) = Q a R2b(q)(Pa — Note that Si will not be a submanifold, because 

is not a regular value of the functions b — m + 1, n. Anyway, the geometric context 
we have developed can be very useful to tackle the problem of the comparison of the two 
methods. 
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Proposition 6.2 If c(t) = {q A {t)) is a solution of the free problem which verifies all the 
constraints, i.e, 

q a (t) = m*(q(t))q a (t) , 1 < a < m, 
then c(t) is a solution of the nonholonomic and vakonomic problems simultaneously. 

dL 

Proof : Let us consider the submanifold 5* = \v n = i*— — |, which is contained in S\. 

1 dq a 

A natural question is whether the vakonomic vector field will be tangent to S, that is, 
X v /~ G TS. ^From equations (0), we have along any integral curve of the vakonomic vector 
field 

d ( JdL\ „ . . A d 2 L ..„ d 2 L 



X vk eTS^-[ Pa - i*— = ^ p a = q A ^— + q 



dt V dq a J ^ dq A dq a * dq a dq a 

On S, we have that 

dL d^ 13 _ dL dL d^ 3 _ dL 
Pa ~ dq^~ Pp ~dq« ~ dq« ~ W~dq^ ~ dq* ' 

Then the above condition can be rewritten as 

dL d f dL 



dq a dt \dq 

that with the third set of equations in (|7|) are precisely the Euler-Lagrange equations. Then, 
we have proved that c(t) is a free solution and satisfies the constraints if and only if 

dT 

is a solution of the vakonomic equations (0). Since the constraints = are automatically 
satisfied for all the points in S we deduce that c(t) is also a solution of the nonholonomic 
problem. \qed\ 



Remark 6.3 As a consequence of Proposition |6.2| we obtain that if g is a riemannian metric 
on Q with kinetic energy L = |<? and we assume that we are given a distribution T> on Q 
which is geodesically invariant with respect to the Levi-Civita connection V 9 , then all the 
nonholonomic solutions can be seen as vakonomic ones. In fact, they all are solutions of the 
free problem. This last result was first stated in || (Theorem 3.2) with additional hypothesis 
on the nature of the metric g and the integrability of V ±9 which are not essential, as we 
have seen. 



Remark 6.4 Let Q : G x Q — ► Q be a free and proper action on Q. Then n : Q — ► Q/G 
is a principal G-bundle. Assume that the lagrangian L : TQ — > K is G-invariant and is 
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subject to equivariant affine constraints, M, such that its linear part T> is the horizontal 
distribution of a principal connection 7 on it : Q — ► Q/G. Then, we have the following 
result, which is an adaptation of Theorem 3.1 in to our geometric description of vakonomic 
and nonholonomic mechanics. 

Proposition 6.5 Assume that the admissibility and compatibility conditions hold. Then, 
the following are equivalent: 

(i) the nonholonomic solution (q A (t),q a (t)) G M verifies the condition gb(q A (t),q a (t),po) = 
for some po, m + 1 < b < n. 

(ii) the curve (q A (t), q a (t),po) G W\ is a vakonomic solution. 



Example 6.6 (Rolling penny [Q) Consider a vertical penny constrained to roll without 
slipping on an horizontal plane and free to rotate about its vertical axis. Let x and y denote 
the position of contact of the disk in the plane. The remaining variables are 9 denoting the 
orientation of a chosen material point P with respect to the vertical and <f> the heading angle 
of the penny. The configuration space is then Q = R 2 x S 1 x S 1 . The la grangian defining 
the dynamical problem may be written as 



1 -(x 2 + y 2 + 9 2 + ^ 2 ) 



2 

and the constraints are given by 

x = 9 cos (j> , 
y = #sin0 . 

For simplicity, we assume that the mass m, the moments of inertia J, J and the radius of 
the penny R are 1. 

We have well defined vector fields 

d ' d ' d \ ' / \ 
X vk = x— + y— + 9— + (p— + -(p [p x sin (p cos (j) - p y cos 2 cf) -29 sin (p) — 

ox ay 39 ocp 2 V / ax 

1 • / \ d 1 • d 
+ -<j> \Px sin 2 (p — p y cos (p sin (p + 29 cos (p ) — + -<p (p x sin (p — p y cos (p) — 7 

2 V ) ay 2 39 

d 

dp x 



d ' / \ d 

—9 (p x sin (p — p y cos (p) — r + (p ( p x sin cp cos (p — p y cos 2 (p — 29 sin 

d(p ^ 

+(p (p x sin 2 (p — p y cos cp sin cp + 29 cos 



dp y 
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on W\ and 



d d ■ d ■ d ■ ■ d ■ ■ d 
X n h = x— + y— + 9— + 0— - 9 sin 0— + 9 cos 0— 
ox oy o9 ocp ox oy 



on M, respectively. Thus, we find that S\ is determined by the vanishing of 

gi = (p x sin - p y cos 0) , 
#2 = -^(^sin0-p y cos0) . 

Si has two connected components 

Ci = {w G Wi | p x sin — p y cos = 0}, 

c 2 = {weWi\e = o, = o}. 

Applying the algorithm, we obtain that Ci = Cu U C12 and (72 = 6*2, where 

On = {weCi\<j> = 0}, 

C12 = {w G C\ I Id = p x cos0 + Py sin0} . 

One step more yields that Cu = Cu and Cu = Cu, so they are also final constraint 
submanifolds. The nonholonomic solutions that fall into Cu are motions of the penny along 
a straight line in the horizontal plane. The nonholonomic solutions in C2 are stationary 
positions. However, any nonholonomic solution can be seen as a vakonomic one contained 
in the final constraint manifold C12, with Lagrange multipliers p x = 2#cos0 and p y = 
28 sin 0. In terms of the extended lagrangian formalism mentioned in Remark [2.4| , we have 
the following Lagrange multipliers 

K = -7. p x = x - p x = -9 COS , 

ox 

^ dL ■ a - a, 

y = ~dy ~ Py = V ~ Py = Sm ' 

which is just the result of Bloch and Crouch Pj. 



Example 6.7 (Constrained particle [ |36|| ) We will discuss here an instructive example 
which has been extensively treated in the literature of nonholonomic mechanics. Consider 
a particle of unit mass moving in space, Q = IR 3 , subject to the constraint 

$ = z — yx = . 

The lagrangian is 

L= 1 -(x 2 + f + z") . 
29 



The constraint is linear and the compatibility and admissibility conditions are satisfied, so 
we have well defined vector fields 

.9 . / d . d \ . d y(p z — z) — yyx ( d d\ . d 

dx \dy dz J dz 1 + y 2 \dx dz J ' dy 

on W\ and 

. d . / d . d \ . d yyx ( d d 

dx \dy dz J dz 1 + y 2 \dx dz 

on M, respectively. Comparing them, we find that Si is determined by the vanishing of 

91 = x(p z - z) , 

92 = -y(p z - z) ■ 

Consequently, Si has two connected components 

Ci = {we Wx\ Pt -z = 0}, 
C 2 = {w E Wi | x = , y = 0} . 

Applying the algorithm, we obtain that Ci = Cn U C12, where 

C n = {weCi\y = 0}, 
C u = {w e Ci I x = 0} . 

On the other hand, C2 = C2, so C2 is a final constraint submanifold. Another step of the 
algorithm yields 

Cu = Cu , C12 = C12 , 
so they both are also final constraint submanifolds. 

Therefore, the nonholonomic solutions that can be seen as vakonomic ones are the ones 
which belong to 



C 



11 



C12 

c 2 



(x t + x , y , z t + z ) , where z = y x 



(x , y t + ?/o, z ) , 
(x , y , z ) , 

that is, stationary or free motions in M. Observe that there are plenty of nonholonomic 
solutions that can not be seen as vakonomic ones. 



Example 6.8 (Ball on a rotating table ||23|| ) Applying the algorithm to this example, 
one can obtain the same result found in [p3| . The configuration space is Q = M 2 x SO (3) with 
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coordinates (x,y,R). We denote the spatial angular velocity by £ G M 3 , where | = AR T . 
The lagrangian is 

^ = ^((e i ) 2 +(a 2 +(a 2 ) + ^(x 2 +y 2 ), 

where I and m are the inertia and mass of the ball, respectively. The constraints are 

x = r£ 2 — VLy , 
y = — r^ 1 + Vtx , 

where r is the radius of the ball and Q is the angular velocity of the table. 
Applying the algorithm, one finds the following final constraint submanifolds 

C fl = {w E Wi | x = y = p x = p y = 0} , 

C f2 = {w e Wi | £ 3 = n} . 



As proved in [23], there are nonholonomic solutions that can not be seen as vakonomic ones. 
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